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We use the Markov Chain Monte Carlo method to investigate a global constraints on generalized 

O^ ' holographic (GH) dark energy with flat and non-flat universe from the current observed data: the 

^— ^ ' Union2 dataset of type supernovae la (SNIa), high-redshift Gamma- Ray Bursts (GRBs), the obser- 

Cn ' vational Hubble data (OHD), the cluster X-ray gas mass fraction, the baryon acoustic oscillation 

jrt ' (BAO), and the cosmic microwave background (CMB) data. The most stringent constraints on GH 

^>. ' model parameter are obtained. In addition, it is found that the equation of state for this generalized 

^vj ' holographic dark energy can cross over the phantom boundary Wdii = — 1. 
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j^ ■ The late accelerating universe p| is often interpreted by introducing a new component dubbed as dark energy 

(DE) with negative pressure in the standard cosmology. And a natural candidate of DE is positive tiny cosmological 

►^ ■ constant, though it suffers from both the fine tuning and cosmic coincidence problems. If DE is not a constant but a 
C_^ I time variable one, the fine tuning and cosmic coincidence problems can be solved. So, lots of dynamical dark energy 
0^ \ models were investigated in the past years [2] ■ Especially, the energy density given by basing the holographic principle 

rn 

O 



I. Introduction 



are studied extensively [3- According to the holographic principle it is required that the total energy for a system 

f^ . with size L should not exceed the mass of a black hole of the same size. The largest L allowed indicates an energy 

density pA = 3c^MpL~^, where c is a numerical constant and Mp is the reduced Planck Mass M^'^ = SirG. Applying 

this principle to cosmology, the UV cut-off is related to the vacuum energy, and IR cut-off is related to the large scale 

^^ ' of the universe such as Hubble horizon, future event horizon, particle horizon, etc. And an accelerated universe can be 

S [ gotten by taking the future event horizon as an IR cut-off, with existing a causality problem. Unfortunately, though 



the Hubble horizon is the most natural cosmological length scale, non-accelerated universe can be obtained [3| when 
this horizon is taken as the IR cut-off. So, how to obtain an accelerated expansion by using the Hubble horizon as 
the IR cut-off is interesting. 

In addition, on the basis of holographic principle Ref. [J| take the Ricci scalar as the IR cut-off and obtain a new 
form of energy density, pj^ — 3c^Mp{H + 2H^ + k/a^) ex R, dubbed as Ricci dark energy. For this model it avoid the 
causality problem and solve the coincidence problem [j|. And in Ref. [5] it is found that the Ricci dark energy has 
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relation with the causal connection scale R~^ = Max{H + 2H^, —H) for a flat universe. Also, it is shown that for these 
two cases only taking R~^ = H + 2H^ as the IR cut-off, the obtained model is consistent with the current cosmic 
observations when the dark energy is looked as an independently conserved component pde + 3H{pde + Pde) — 
[5|. And as indicated in Ref. [6|, H^ or H alone can not provide an late accelerated universe that is consistent 
with the current cosmic observations. So, the generalized holographic (or Ricci) dark energy model, i.e. a form of 
their combination is investigated in Ref. [6[. In this paper we applying the current observed data to constrain the 
generalized holographic (GH ) dark energy model by using the Markov Chain Monte Carlo (MCMC) method. 

II. Basic equations for generalized holographic dark energy 

In a Friedmann-Robertson- Walker universe, when the Hubble horizon and Ricci scalar are taken as the IR cut-off, 
the holographic dark energy and Ricci dark energy are written as, ph — 3c^MpH^ and pn = ic^MpR, respectively. 
And in order to compare the holographic and the Ricci dark energy, and to obtain an accelerated universe by using 
the Hubble horizon as the IR cut-off, in Ref. j6(] a generalized version of holographic dark energy are constructed as, 

PGH = icHlfj{^)H\ (1) 

where /(x) is a function of the dimensionless variable x = R/H^, and it is interesting to write the function as |6|, 

where e is parameter. For the generalized form of energy density, when e = or e = 1, it becomes holographic or 
Ricci dark energy density, respectively. Thus for this generalized model, the dark energy density is expressed as 

PGH = 3cHl^[l - e{l - ^)]H^ 

o Or . H -\- zH , ^ o 

= 3cX'[l + e-e(l + .)l^]i/2. (3) 

And its dimensionless dark energy density is described, 

o = Peg 



9r-, /, dh\H 

c2[l + e-e(l + z)— — ]. (4) 

az 



For the generalized holographic dark energy model, the corresponding Friedmann equation can be written as, 

rr2 _ rr2X^Gm{l + zf + nr{l + zY + ^k{l + zf) 2(f]om + f^r + ^fc) w-, , .2--^ + ^^ ,.. 

^ — ^ov TT"; — 57 — ?v; ^ \^ T"; — 57 — ?r; — JU + ^J " ' N' v^) 



where r^om, ^r and ^k respectively denotes the current value of dimensionless matter, photon and curvature density, 
here S^om include baryon matter f2f, and cold dark matter U.^, f^om = ^6 + ^c- Furthermore, for the deceleration 
parameter q[z) and the geometrical diagnostic quantity Om{z) 7], they can be expressed by the Hubble parameter 
as, 

" , ,1 dH 

-i + (i + ^)T7-jr> (6) 



am ^ ' H dz' 



3 

Om{z) = \\ , x^l + z = -, E{z) = -^. 7 

And the equation of state (EOS) of generalized holographic dark energy wgh can be derived as, 

(1 + z) 1 dpGH ,gx 

WGH = -l + — 5 1 — , (8) 

3 pgh dz 

according to the conservation equation with no interactions between two dark components pgh+'^H{1+wgh)pgh = 0. 



III. The current observed data and cosmological constraint methods 

In this part we introduce the cosmological constraint methods and the current observed data used in this paper. 
Concretely, it includes 557 Union2 dataset of type supernovae la (SNIa) [8], 59 high-redshift Gamma- Ray Bursts 



(GRBs) data [9(, observational Hubble data (OHD) 10], X-ray gas mass fraction in cluster [ll|, baryon acoustic 
oscillation (BAO) 12*, and cosmic microwave background (CMB) data [13| . 



A. Type la supernovae 

For SNIa observations, we use the SNIa Union2 dataset that includes 557 SNIa [8]. Following [ij, ll5[, one can 
obtain the corresponding constraints by fitting the distance modulus p{z), 

Hthiz) ^5\ogj^o[DL{z)]+ fiQ. (9) 

In this expression Dl{z) — iJo(ii(z)/c is the Hubble-free luminosity distance, with Hq being the Hubble constant 
described by the re-normalized quantity h as Hq = lOOh km s^^Mpc^^, and 

Cjl + Z) . r—- [^ dz' 

dL[z) = — j==r s^nn[^/ \Uk\ 



TT-l 

Mo - 5Zogio(-r^) + 25 = 42.38 - Uog^h, 
Mpc 



where sinnn{^J\Q.k\x) respectively denotes sin(-\/|fTfe|x), y^fTfc|x, and sinh(-\/|STfe|a;) for ilfc < 0, J7fc = and Q.k > 0. 
Additionally, the observed distance moduli p.ohs{zi) of SNIa at z^ is 

tJ^obs[zi) = mobs[zi) - M, (10) 

where m and M are apparent magnitude and absolute magnitude of SNIa. 

For using SNIa data, theoretical model parameters ps can be determined by a likelihood analysis, based on the 
calculation of 

SNIa ^^ 

sp {^^ogiQ[DL{Ps,Zi)] - mobsjzi) + M'}^ 

SNIa ^i 

where M' = /.iq + M is a nuisance parameter which includes the absolute magnitude and the parameter h. The 
nuisance parameter M' can be marginalized over analytically |lq] as 



cxp 

-CXD 



+ 00 



-\x^{ps,M') 



resulting to 



^''^'^^'"{tJ- (1^' 



with 

1 2 



SNIa * 

^ ^? ' 

SNIa * 
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Relation (jlip has a minimum at the nuisance parameter value M' — B/C, which contains information of the values 
of h and M. Considering that the expression 

xlNiaiPs)=A-iByCl (13) 

is only different from Eq. (fT2|) with a constant term ln(C/27r), it is often used in the likelihood analysis [14Lll6l|. 

B. High-redshift Gamma- Ray Bursts data 

The GRBs data can be observed at higher redshift than SNIa. The currently observed reshift range of GRBs is at 
0-^ ^ z < 9. Therefore, the GRBs data can be viewed as an excellent complement to SNIa data and would provide 
more information at high redshift. When several empirical relations of the GRBs are proposed, these indicators have 
motivated the authors make use of the GRBs as cosmological standard candles at high redshift. However, the fact that 
there are not sufhcient low reshift GRBs leads that the calibration of GRB relations is dependent on the cosmological 
model, namely, the circularity problem. One of methods to solve the circularity problem is the calibration of GRB 
relations are performed by the use of a sample of SNIa at low redshift in the cosmology-independent way [17|. Here, 
the GRBs data we used consists of 59 GRB samples with a redshift range of 1.4 < z < 9 obtained in 9]. These 59 
GRBs are calibrated by utilizing the newly released 557 Uion2 SNIa and the isotropic energy-peak spectral energy 
{Eiso- Ep,i) relation (i.e. Amati relation) 18]. 

The Xgrbs takes the same form as XsNia 

2 . N Mr-'[fJ-obs{zt- fJ.th{zt;Ps,fJ.o)? ,,.. 

XGRBsiPs,tJ-a) = 2^ -^ • (14) 

1=1 * 

The same method are used to deal with the nuisance parameter hq as shown in the description of XsNia above. 

C. Observational Hubble data 

The observational Hubble data [19| are given by basing the differential ages of the galaxies. In 20], Jimenez et al. 
obtain an independent estimate for Hubble parameter, and use it to constrain the cosmological models. The Hubble 
parameter as a function of redshift z can be written in the form of 

H{z)^--^^. (15) 

1 + z at 



So, once dz/dt is known, H{z) is obtained directly. By using the differential ages of passively-evolving galaxies, Refs. 
tld I21I, 122] obtain twelve values of H{z) at different redshift (redshift interval < z < 1.8), as listed in Table HI In 



z 


0.1 0.17 0.27 0.4 0.48 0.88 0.9 1.30 1.43 1.53 1.75 


H{z) (kms-^Mpc-i) 


74.2 69 83 77 95 97 90 117 168 177 140 202 


la uncertainty 


±3.6 ±12 ±8 ±14 ±17 ±60 ±40 ±23 ±17 ±18 ±14 ±40 



TABLE I: The observational H{z) data [2l|,|2 

addition, in [23| the authors take the BAO scale as a standard ruler in the radial direction, and obtain three more 
additional data: H{z = 0.24) = 79.69 ± 2.32, H{z = 0.34) = 83.8 ± 2.96, and H{z = 0.43) = 86.45 ± 3.27 



The values of model parameters can be determined according to the observational Hubble data by minimizing 24 1 

Xohd{Ho,Ps) = 2^ JT^; , (16) 

where Hth is the predicted value of the Hubble parameter. Hobs is the observed value, cr(zi) is the standard deviation 
measurement uncertainty, and the summation is over the 15 observational Hubble data points at redshifts Zi. 



D. The X-ray gas mass fraction 

The X-ray gas mass fraction, fgas, is defined as the ratio of the X-ray gas mass to the total mass of a cluster, 
which is approximately independent on the redshift for the hot (fcT > 5keV), dynamically relaxed clusters at the 
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, the ACDM model is much favored and is chosen as 



radii larger than the innermost core r25oo- As investigated in 

the referenced cosmology. The model fitted to the referenced ACDM data is presented as [Hj 



fACDM. N _ KA-fb{z) ( Ob 



/^a"."^'^(^) 



Da{z) 



(17) 



where D^'^^^^ {z) and Da{z) denote respectively the proper angular diameter distance in the ACDM cosmology and 
the current constraint model. A is the angular correction factor, which is caused by the change in angle for the current 
test model 6*2500 in comparison with that of the reference cosmology ^'^500'*^: 



92500 ; \[H{z)Da{z)]^^d^^ 

here, the index 77 is the slope of the fgas{'r/r25oo) data within the radius r2500 7 with the best-fit average value 
77 = 0.214 ± 0.022 11]. And the proper (not comoving) angular diameter distance is given by 



"'''' = iiT^K^^-i^y. 7?m'- '''' 



It is clear that this quantity is related with d^i^z) by 

dL{z) 



Da{z) 



{l + zf 

In equation (|17p. the parameter 7 denotes permissible departures from the assumption of hydrostatic equilibrium, due 
to non-thermal pressure support; the bias factor h{z) — 60 (1 + Ubz) accounts for uncertainties in the cluster depletion 



factor; s{z) = sq{1 + agz) accounts for uncertainties of the baryonic mass fraction in stars and a Gaussian prior for 



So is employed, with so — (0.16 ± 0.05)h'^Q Uj; the factor K is used to describe the combined effects of the residual 
uncertainties, such as the instrumental calibration and certain X-ray modelling issues, and a Gaussian prior for the 
'calibration' factor is considered by ^ = 1.0 ± 0.1 [11]. 

Following the method in Ref. Uj, |25| and adopting the updated 42 observational fgas data in Ref. [lj|, the values 



of model parameters for the X-ray gas mass fraction analysis are determined by minimizing, 

2 _^ [f^a's''^'i^^)-f9as{zi)? , (sp-O.lB)^ , (X-1.0)^ , (ry - 0.214)^ 
XcBF-2^ ^2 ( ) + 0.00162 + 0.012 + 0.0222 ' ^^^) 

i Jgas ^ 

where af ^^(^i) is the statistical uncertainties (Table 3 of 11]). As pointed out in [ll| . the acquiescent systematic 
uncertainties have been considered according to the parameters i.e. 77, b{z), s{z) and K. 



E. Baryon acoustic oscillation 

The baryon acoustic oscillations are detected in the clustering of the combined 2dFGRS and SDSS main galaxy 
samples, which measure the distance-redshift relation at zbao — 0.2 and zbao — 0.35. The observed scale of the 
BAG calculated from these s amp les, are analyzed using estimates of the correlated errors to constrain the form of the 
distance measure I?y(z) 12ll26j 



Dviz) = [(1 + zfDliz)-^]'^' = Ho[-^i r ^;^^f]-. (21) 

H[z;ps) E{z;ps) Jq E{z ;ps) 

In this expression E{z;ps) — H{z;ps)/Ho. The peak positions of the BAG depend on the ratio of Dv{z) to the sound 
horizon size at the drag epoch (where baryons were released from photons) Zd, which can be obtained by using a 
fitting formula 

with 

61 = 0.313{nomh^r"'^^^[l + Omiiriomhy-^''*], (23) 

62 = 0.238(f]om/i^)°-^^^ (24) 

In this paper, we use the data of rs{zd)/Dv{z) extracted from the Sloan Digitial Sky Survey (SDSS) and the Two 
Degree Field Galaxy Redshift Survey (2dFGRS) !2d|, which are listed in Table HTl where rs{z) is the comoving sound 
horizon size 

' " dt r Csda r°° , Cs 



rs{z) =c — = c -T77 = c/ dz 







a 



a^H 7, H{z) 



c /■i/(i+^) da 

' (25) 



V3 7o a^H{a)^yi + {'3nb/{'in^)a) 

where Cg is the sound speed of the photon— baryon fluid 



c- = 3+^x^.3 + lx(§^)a. (26) 

3 P-y\Z) 3 SZ-), 



z 


r,{za)IDv{z) 


0.2 


0.1905 ±0.0061 


0.35 


0.1097 ±0.0036 



TABLE II: The observational rs{zi)IDv{z) data 



Ea^ 



Using the data of BAO in Table |TT] and the inverse covariance matrix V ^ in 12[ : 

/ 30124.1 -17226.9 \ 
V-^=\ , (27) 

\^ -17226.9 86976.6 j 

tlie Xbao(Ps) is gi'^en as 

x\ao{Vs) = X'V-'X, (28) 

where X is a column vector formed from the values of theory minus the corresponding observational data, with 



rs(zd) 



0.1905 



Dv/(0.35) 

and AT* denotes its transpose. 

F. Cosmic microwave background 

The CMB shift parameter R is provided by 



X = ^-/of - - - - I , (29) 

' Ii(£iiL _ 0.1097 ' 



Hgdz 



R = JnornH^il + z,)Da{z.)/c = V^W / -nr-, r, (30) 

here, the redshift z.^, (the decoupling epoch of photons) is obtained using the fitting function 

z* = 1048 [1 ± 0.00124(fi6;i2)-o.738j ^^ _,_ g^(yiamhy''] , 

where the functions gi and (?2 read 

31 = 0.0783(f}b/i2)-"-238(l + 39.5(r!fc/i2)0.763)-i^ 

32 = 0.560 (1± 21. l(f}b/(,2)i-8i)"\ 

In addition, the acoustic scale is related to a distance ratio, I?yi(z)/rs(z), and at decoupling epoch it is defined as 

lA^il + z.f-^. (31) 

where Eq. (j31l) arises a factor 1 + z,, because Da{z) is the proper angular diameter distance, whereas r<,(z») is the 
comoving sound horizon. Using the data of lA,R,Zf, in [13j and their covariance matrix of [Ia{z^,),R{z^),z^] (please 
see table Hill and IIV[) . we can calculate the likelihood L as Xcmb = — 21nL: 

xImb = AddCov-\d,,d,)[Ad,]% (32) 

where Adi — di — df""*""" is a row vector, and di — [Ia, R,z^). 
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FIG. 1: The 2-D contours with Icr, 2a confidence levels and 1-D distribution of model parameters in the non-flat GH model. 
Solid lines are mean likelihoods of samples, and dotted lines are marginalized probabilities for ID distribution. 



7-year maximum likelihood error, a 



Ia{z,) 


302.09 


0.76 


R{z,) 


1.725 


0.018 


z. 


1091.3 


0.91 



TABLE III: The values of Uiz-^), R{z,), and 2. from 7-year WMAP data. 



IV. Observed constraints on generalized holographic DE model by using MCMC method 



Next we apply the Markov Chain Monte Carlo method to investigate a global constraint on above generalized 
holographic dark energy model. The MCMC source code can be found in the CosmoMC package [27| and the 
modified CosmoMC package [uJ, l28|, l29| (this package is about the constraint code of X-ray cluster gas mass fraction) . 
To get the converged results, in MCMC calculation we test the convergence of the chains by taking i? — 1 to be less 



Ia{z,) R(z,) z, 

Ia{z,) 2.305 29.698 -1.333 
R{z.^) 6825.270 -113.180 

2* 3.414 

TABLE IV: The inverse covariance matrix of /a(z*), R{z-t), and z, from 7-year WMAP data. 



0.022 0.025 
0.14 

0.12 



0.022 0.025 0.12 0.14 



0.55 vi:>/ 0.55 \iri> 



0.022 0.025 0.12 0.14 0.55 0.6 

2 



1.5 



^ 0.75 

o 

a° 0.7 



1.5 



1.5 



0.022 0.025 0.12 0.14 0.55 0.6 1.5 2 

0.75 r 



0.75 
0.7 



0.75 
0.7 



0.022 0.025 0.12 0.14 0.55 0.6 1.5 2 



S 13.5 
S' 13 



14 

13.5 

13 



14 

13.5 

13 



14 

13.5 

13 



0.7 0.75 



14 

13.5 

13 



0.022 0.025 0.12 0.14 0.55 0.6 1.5 2 0.7 0.75 13 13.5 14 

.° °'^ °"^ °-'^ °'^ °"^ °"^ 
0.25 I I 0.25 I' — I 0.25 I ^ — <-J 0.25 I ^-^ I 0.251 3iJ 0.25 L^^^ 



0.022 0.025 0.12 0.14 0.55 0.6 1.5 2 0.7 0.75 13 13.5 14 0.25 0.3 



75 



70 



75 



70 



75 



75 



70 



75 



75 



70 



75 



70 



0.022 0.025 0.12 0.14 0.55 0.6 1.5 2 0.7 0.75 13 13.5 14 0.25 0.3 



n h^ 



Age/Gyr n^^ 



70 75 
H„ 



FIG. 2: The 2-D contours with la, 2a confidence levels and 1-D distribution of model parameters in the fiat GH model. Solid 
lines are mean likelihoods of samples, and dotted lines are marginalized probabilities for ID distribution. 



than 0.03. The total y^ is expressed as, 

2/\_2 |2 |2 |2 |2 |2 

XtotalKPsj — XsNIa + XgRBs + XoHD + XcBF + XbAO + XcMBi 

with the parameter vector reading 

Ps = {nbh^,nch^,nk,e,c}. 



(33) 



(34) 



Here the expression of x^ for each observation corresponds to Egs. p^ . HH), (HH), (PO)) . ([^5]) and (|32p . Based on the 
basic cosmological parameters Ps we can also obtain the derived parameters flom = ^b + ^c, ^ogh = 1^ ^Om ^^k, and 
the Hubble constant Hq — lOO/i km-s ^^-Mpc"^. Using the currently observed data with the Xtotai ™ ^'i- P^)) - Figs. 
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[T] and [2] plot the 2-D contours with la, 2cr confidence levels and 1-D distribution of model parameters in the flat and 
non-flat generalized holographic dark energy model. Solid lines are mean likelihoods of samples, and dotted lines are 
marginalized probabilities for ID distribution. Table IVl lists the MCMC calculation results for the constraint on model 
parameters. It includes the means, standard deviations with the marginalized limits for the model parameters, and 
the values for the best-fit sample, and projections of the n-dimensional la and 2a confidence regions. The n-D limits 
give some idea of the range of the posterior, and are much more conservative than the marginalized limits 27[. From 
the table |V] it can be seen that for the non-flat universe, the best fit results are given as flk = ~0.0047lg'Qj!)gglQ'gj^20' 
c = 0.576to mt-Vall^ « = '^■^'^^-oitl-OA^ ^Om = 0.280l°;°^^t°°^^ (it has a smaller value of t^om relative to the case 
of constraints on the Ricci dark energy model 30|, |3l|), with Xmin — 619.314. And for this case, it predicts the age 



of universe tage — 13.71lto[slgt.o'^^si^'y^)- Furthermore, comparing the Ref. [6| one can see that for the generalized 
holographic dark energy the more stringent constraint on model parameters at 2a confidence level are given in this 
paper by using the more observational data, and it tends to have a smaller value of dimensionless matter density fiom 
and a bigger value of model parameter e. 





Non-flat 


Non-flat 


Flat 


Flat 


Parameters 


Best fit values 


Means 


Best fit values 


Means 


Qth^ 


n n9'?'?+00023+0.0027 
"•'-'^■J'J-0. 0016-0. 0016 


n n9'^K+0'0006+0.0013 
U.UZJD_o. 0006-0. 0012 


n n9'?R+00018+0.0022 
U.UZOD_o 0017-0.0022 


n noqe+o. 0006+0. 0012 
u.uzjD_o. 0007-0. 0012 


n^h" 


n 1 1 cn+00220+0.0290 
U.iiOU_o. 0160-0. 0160 


n 1 1 00+0. 0067+0. 0133 
U.iiOO-0. 0065-0. 0120 


n 1 1 7S+0. 0195+0. 0241 
U.ii/ 8-0.0105-0.0121 


n -10-17+0.0055+0.0126 
U-J^J'-0. 0056-0. 0097 


Vtk 


n nnzL7+0'0132+0.0159 
u.uu-il -0.0089-0. 0120 


n nn9Q+0'004o+o.oo77 

U.UUZy_o 0040-0.0168 
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c 


n cyc+0. 034+0. 037 
U.OrD_Q o3g_o 053 


n c;7/|+0. 013+0. 024 
U.O(y:_o 012-0.026 


n CSR+0. 023+0. 027 
U.080_o 043-0.052 


n c;7c;+0. 012+0. 023 
U.OrO_o 013-0.025 


e 


1 O/lQ+0. 347+0. 461 
r.o'±y_o. 380-0. 442 


1 SI c;+0. 126+0. 262 
i.OiO_o i29_o.237 


1 S/iq+0-392+0.461 
4.81O-0. 347-0. 429 


1 snS+0-130+0.262 
±.OUO_o 129-0.255 


^Om 


n 9on+0-036+0.047 
U.ZOU_o 032-0.040 


n 9S1 +0.013+0.026 
U.ZOi_o 013-0.024 


n 97Q+0. 034+0. 051 
U.Z/y_o 035-0.037 


n 9S'?+0.013+0.029 
U.ZOO_o 013-0.024 


^OGH 


n 79c;+0. 035+0. 041 
U. '^>J_o. 044-0. 055 


n 790+0.013+0.025 
U.(ZZ_o 014-0.028 


r, 79-1+0.035+0.037 
U. 'ZJ-0. 034-0. 051 


n 717+0.013+0.024 
U-' J '-0. 013-0. 029 


Ho 


7n Ofii +4.710+5.143 


71 1 SI +1.313+2.627 
/i. 181-1 274-2.490 


71 170+3.514+3.939 
ri. 1^-2. 369-3. 347 


71 RKfi+l. 102+2. 056 
^i.DOD-i 095-2.151 


tageiGyv) 


io 71 1+0.709+0.924 
10. 1 J-J--o.859-0.978 


1Q C/lo+0.263+0.524 
ij.04y_o. 267-0. 515 


-10 ..<(;r,+0. 565+0. 584 
iOi.4DZ-o. 589-0. 747 


1 ■? /ll 0+0. 198+0. 381 
10.^110-0.200-0.402 



TABLE V: For the fiat and non-fiat universe, the best fit model parameters with their limits from the extremal values of the 
n-dimensional distribution (recommended); and the means with the marginalized limits for the model parameters, from MCMC 
calculation by using SNIa Union2, GRBs, OHD, CBF, BAG, and CMB data. 





zt 


90 


Omo 


WOGH 


Non-fiat 


0-706l'o.lll 


-0.639lHf7 


O.241IHI7 


1 nc;i+oo48 

— i.UOi_o 048 


Flat 


0.705t°:«^« 


-O.SQStHl^ 


0.2681H1^ 


1 ni c;+oo45 

i.UiO_o.o45 



TABLE VL The best fit values of transition redshift, current values of deceleration parameter, Om parameter, and EOS of 
generalized holographic dark energy with their confidence levels for flat and non-flat universe. 



In addition, according to the calculation of the covariance matrix and the best fit values of model parameters, 
the best fit evolutions of deceleration parameter q{z), geometrical quantity Om{z) and EOS of dark energy wgh{z) 
with their confidence level (shadow region) are plotted in Figs. [3]and[4l From the figures we can see that a current 
accelerated universe is obtained, and the equation of state for this generalized model can cross over the boundary 
of cosmological constant Wf^{z) — —1. And for this generalized dark energy model, the predicted values of some 
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K 


ho 


Oib 


So 


Qs 


Non-flat 


0.212 


1.081 


0.998 


0.732 


-0.092 


0.174 


-0.055 


Flat 


0.208 


1.025 


0.958 


0.783 


-0.086 


0.156 


0.020 



TABLE VII: The best fit values of parameters in fgas analysis method for flat and non-flat universe. 
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FIG. 3: The evolution of q{z), Om[z) and Wdeiz) for non-flat generalized holographic model. 



cosmological parameters with flat and non-flat universe according to above combined constraints are listed in table 
IVII From this table, it can be found that the current values of the deceleration parameter and the EOS of GH model 
are, qo = -0.639+[J;J]^^, w^gh = -l.OSltaoS for non-flat universe, and go = -0.598+[J;J]^^, woGff = -IMbt^^lf^ ^r 
flat universe. At last as an appendant, in table IVIII we also show the best fit values of several parameters in fgas 

analysis method. 

. . - 1 . . 

By the way, in appendix we also list the constraint results on another generalized model in Ref. i6| , i.e. generalized 

Ricci DE by using the MCMC method and above observed data. 



Conclusions 



In summary, for interpreting the accelerating universe and solving the coincidence problems of cosmological constant, 
the holographic dark energy models are extensively studied from the different points of view. In holographic cosmology, 
considering that taking the natural Hubble horizon as the IR cut-off to obtain an accelerated universe is interesting, 
Ref. y presents a new generalized holographic dark energy model. In physics, this generalized model investigate a 
new idea to interpret the accelerating universe by using the holographic principle with including the Hubble horizon as 
an IR cut-off. In addition, the holographic and Ricci dark energy can be compared in the generalized model according 
to the new introduced parameter e. In this paper, the flat and the non-flat generalized holographic dark energy are 
constrained according to the current observed data. The stringent constraints on model parameters are given from the 
MCMC calculation. Considering the cosmic constraint on the parameter e, it is obtained that the cosmic data favor 
a generalized dark energy model which is more Ricci-like, since one has the relation pcff = ^Pr. + (1 ^ ()ph and the 
best fit value of parameter e = 1.849 for a non-flat universe constrained from the observational data. And according 
to the constraint results, it is shown that relative to the Ricci dark energy model {^om — 0.300lg Q37]']gQ42 3l|), 
it has a smaller value of the dimensionless matter density J7om = 0-280lg 032^0040 f*-"^ ^^'^ non-flat universe, which 
result is more consistent with the current observations and cosmological constant model ;32]. In addition, based on 
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^0.29 

d 0.27 

0.25 
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FIG. 4: The evolution of q{z), Om{z) and Wde{z) for flat generalized holographic model. 

the calculation of covariance matrix the best fit evolutions of cosmological quantities such as deceleration parameter, 
Om parameter and EOS of generalized holographic dark energy with their confidence region are discussed. It is found 
that the EOS for this dark energy model can cross over the boundary of cosmological constant {wj\ = — 1). And the 
values of transition redshift, current deceleration parameter, EOS of GH dark energy are obtained, respectively. It 
can be seen that for the flat universe the best fit value of wogh = ^1-O15too45 is near to the cosmological constant 
model. 
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Appendix A: cosmological combined constraints on generalized Ricci DE model by using MCMC method 

Considering Ref. 6] , another extended form dubbed as generalized Ricci dark energy is expressed. 



ij2 



paH = 3c'M^[l-ij{l-—)]R, 



R 



(Al) 



where 77 is a parameter. It is easy to see when 77 = 1 or 77 = 0, this generalized form reduces to Ricci or holographic 
dark energy, respectively. And the Friedmann equation is described for this generalized model as, 

rj2_rj2u^ 2(»o„ + ». + »fc) ^ ^^ , ^^^—+-An^ , 2{no,^{l + zf + ^rjl + zf + ^kjl + zf) , 

^ -^0^^^ 2-c2(l + r;) ^^^ + ^^ + 2-c2(l+ry) ^^ ^^^^ 

From above equations one can see that two generalized dark energy models are equivalent when e — 1 — -q. Figs. [S] 

shows the \D distributions of model parameters. And the MCMC calculation results for the non-flat universe are, 
O, - -0 nnn8+°™89+°-°^28 r, ,-0^+0.024+0.029 „ _ _« oc-c- +0.374+0.470 - q _ n ofin+ooa^+o.oso r ■- - ■ 

\Lk — U.UUU»_Q Qj27~0.0143' "^ " ^'^"^-0. 045-0. 049' ^ — "■°00_o.366-0.448' ^^'^ "Om — U.Z»U_q q32_o.036 ^^^ ^^^ "^^^ 

fit values. 
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